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Abstract. As several large scale interferometers are beginning to take data at sen- 
sitivities where astrophysical sources are predicted, the direct detection of gravita- 
tional waves may well be imminent. This would (finally) open the long anticipated 
gravitational-wave window to our Universe, and should lead to a much improved 
understanding of the most violent processes imaginable; the formation of black holes 
and neutron stars following core collapse supernovae and the merger of compact 
objects at the end of binary inspiral. Over the next decade we can hope to learn 
much about the extreme physics associated with, in particular, neutron stars. 

This contribution is divided in two parts. The first part provides a text-book 
level introduction to gravitational radiation. The key concepts required for a dis- 
cussion of gravitational-wave physics are introduced. In particular, the quadrupole 
formula is applied to the anticipated "bread-and-butter" source for detectors like 
LIGO, GEO600, EGO and TAMA300: inspiralling compact binaries. The second 
part provides a brief review of high frequency gravitational waves. In the frequency 
range above (say) 100 Hz, gravitational collapse, rotational instabilities and oscil- 
lations of the remnant compact objects are potentially important sources of grav- 
itational waves. Significant and unique information concerning the various stages 
of collapse, the evolution of protoneutron stars and the details of the supranu- 
clear equation of state of such objects can be drawn from careful study of the 
gravitational-wave signal. As the amount of exciting physics one may be able to 
study via the detections of gravitational waves from these sources is truly inspir- 
ing, there is strong motivation for the development of future generations of ground 
based detectors sensitive in the range from hundreds of Hz to several kHz. 



1 Introduction 



One of the central predictions of Einsteins' general theory of relativity is 
that gravitational waves will be generated as masses are accelerated. De- 
spite decades of effort these ripples in spacetime have still not been ob- 
served directly. Yet we have strong indirect evidence for their existence from 
the excellent agreement between the observed inspiral rate of the binary 
pulsar PSR1913+16 and the theoretical prediction (better than 1% in the 
phase evolution). This provides confidence in the theory and suggests that 
"gravitational-wave astronomy" should be viewed as a serious proposition. 
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Provided that i) detectors with the required sensitivity can be constructed, 
and ii) the significant data analysis challenge can be dealt with, this new win- 
dow to the Universe promises to bring unprecedented insights into the most 
violent events imaginable; supernova explosions, binary mergers and the big 
bang itself. A key reason for this expectation follows from a comparison be- 
tween gravitational and electromagnetic waves: 

While electromagnetic waves are radiated by individual particles, gravi- 
tational waves are due to non-spherical bulk motion of matter. In essence, 
this means that the information carried by electromagnetic waves is 
stochastic in nature, while the gravitational waves provide insights into 
coherent mass currents. 

The electromagnetic waves will have been scattered many times. In con- 
trast, gravitational waves interact weakly with matter and arrive at the 
Earth in pristine condition. This means that gravitational waves can be 
used to probe regions of space that are opaque to electromagnetic waves, 
It is, of course, a blessing in disguise since the weak interaction with 
matter also makes the gravitational waves fiendishly hard to detect. 
Standard astronomy is based on deep imaging of small fields of view, while 
gravitational-wave detectors cover virtually the entire sky. 
Electromagnetic radiation has a wavelength smaller than the size of the 
emitter, while the wavelength of a gravitational wave is usually larger 
than the size of the source. This means that we cannot use gravitational- 
wave data to create an image of the source. In fact, gravitational- wave 
observations are more like audio than visual. 

Morale: Gravitational waves carry information which would be very difficult 
to glean by other means. By analysing gravitational- wave data we can expect 
to learn a lot about the extreme physics governing compact objects. This 
should lead to answers to many outstanding questions in astrophysics, for 
example, 

What is the black-hole population of the Universe? Observations of the 
nonlinear spacetime dynamics associated with binary merger, as well as 
the quasinormal mode (QNM) ringing which is likely to dominate the 
radiation at late times, should provide direct proof of the presence of a 
black hole, as well as a measure of it's mass and spin. 
Arc astrophysical black holes, indeed, described by the Kerr metric? By 
studying the inspiral of a low-mass object into a supramassive black hole 
(as anticipated at the centre of most galaxies) we can hope to construct 
a detailed map of the exterior black-hole spacetime. 

What is the supranuclear neutron star equation of state? This is a 
very difficult question, which potentially involves an understanding of 
the role of exotic phases of matter, like deconfined quarks, superfluid- 
ity/conductivity, extreme magnetic fields etcetera. At an immediate level, 
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gravitational-wave observations of the various oscillation modes of a com- 
pact star may be used to solve the "inverse problem" for parameters like 
the mass and the radius of the star. At a more subtle level, the spectrum of 
the star depends on the internal physics. For example, strong composition 
gradients which affect the gravity g-modes significantly, and superfluidity 
which leads to the existence of new classes of oscillations associated with 
the more or less decoupled additional degrees of freedom of such a system. 
Finally, the internal physics also plays a crucial role in determining the 
extent to which various rotational instabilities will be able to grow to an 
interesting amplitude. In the case of the Coriolis driven r-modes, their 
instability may be severely suppressed by the presence of hyperons in the 
stars core. A key role is also played by the internal magnetic fields. 

At the time of writing, the new generation of intcrferomctric gravitational- 
wave detectors (in particular LIGO and GEO600) is already collecting data at 
a sensitivity at least one order of magnitude better than that of the operating 
resonant detectors. In the first instance, the broadband detectors will be 
sensitive in a range of frequencies between 50 and a few hundred Hz. This 
frequency window is of great interest since an inspiraling compact binary will 
move through it during the last few minutes before merger. Such sources are 
the natural "bread and butter" source for the detectors. The next generation 
of interferometers will broaden the bandwidth somewhat but will still not be 
very sensitive to frequencies above 500-600 Hz, unless they are operated in 
a narrow-band configuration [1, 2]. There are also interesting suggestions for 
wide band resonant detectors in the kHz band [3]. The move towards higher 
frequencies is driven by the wealth of exciting sources that radiate in the 
range from a few hundred Hz up to several kHz. 

Our contribution to these proceedings is divided in two parts. The first 
part describes gravitational- wave physics at an introductory level. The second 
part provides a brief review of the main sources that radiate in the frequency 
band above a few hundred Hz. We believe that these sources are the natural 
targets for a third generation of ground based detectors. As we will discuss, 
there are a variety of sources associated with very interesting physics in this 
high-frequency window. These sources clearly deserve special attention, and 
if either resonant or narrow-band interferometers can achieve the required 
sensitivity, a plethora of unique information can be gathered. 

2 Einstein's elusive waves 

The aim of the first part of our contribution is to provide a condensed text- 
book level introduction to gravitational waves. Although in no sense complete 
this description should prepare the reader for the discussion of high-frequency 
sources which follows. 
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2.1 The nature of the waves 

The first aspect of gravitational waves that we need to appreciate is their tidal 
nature. This is important because it implies that they can only be measured 
through the relative motion of bodies. That this should be the case is easy 
to understand. In general relativity we can always construct a local inertial 
frame associated with a given observer. In this local frame, spacetime will 
by construction be flat which means that we cannot hope to observe the 
local deformations which would correspond to a gravitational wave. Consider 



t 




Fig. 1. An illustration of the two test particles discussed in the text. If a gravita- 
tional wave passes through, particle A will observe particle B oscillating back and 
forth. 



two test particles, A and B, that are initially at rest. Assume that they are 
separated by a purely spatial vector £ J (hereafter latin indices are spatial 
and run from 1-3, while greek indices are spacetime and run from 0-4), and 
use the local inertial frame in which particle A remains at the origin for the 
calculation. In this case the equation of geodesic deviation can be written 



dt 2 - K ow 



R'oko^ (1) 



where R^s-y is the Riemann tensor. Here it represents the curvature induced 
by the gravitational wave. Letting £■? = x^ + Sx? , with Sx^ the small displace- 
ment away from the original position, we get 

Q 2 l) X 3 

dt 2 ~ ~ R \ko x o = -RjQkox a (2) 
Now it is natural to define the gravitational wave-field hj k through 
_ 1 d 2 h jk d 2 6x 3 1 d 2 h jk „ fc 



which integrates to 
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5x 3 = 2 h i kX o or h ~ — ( 4 ) 

where h is the dimensionless gravitational-wave strain. From this exercise we 
learn that, in order to detect gravitational waves we need to monitor (with 
extreme precision) the relative motion of test masses. Let us now assume that 




the waves propagate in the ^-direction, i.e. that we have hjk — hjk(t — z). 
Then one can show that we have only two independent components; 

^+ = ^xx = ~hyy h x — h^y — hyj (5) 

What effect does h+ have on matter? Consider a particle initially located at 
(xo, yo ) an d let 7i x = to find that 

Sx = ^h + x (6) 

Sy = --^h+yo (7) 

That is, if h + is oscillatory (a wave!) then an object will first experience a 
stretch in the x-direction accompanied by a squeeze in the t/-direction. One 
half-cycle later, the squeeze is in the x-direction and the stretch in the y- 
direction. It is straightforward to show that the effect of h x is the same, but 
rotated by 45 degrees. This is illustrated in Fig. 2. A general wave will be a 
linear combination of the two polarisations. 

But wait a second! We are discussing the effect of gravitational waves 
without actually having proven that Einstein's theory predicts their existence. 
To remedy this, consider small perturbations away from a flat spacetime. 
That is, use g a p — r) a p + h a p to get the (linearised) Ricci tensor: 
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In vacuum, Einstein's equations are equivalent to requiring that 

Rp. v = (9) 
Now consider what happens if we make a coordinate transformation. For 

/ -» x a + C (10) 

we get 

fo[iis * h^if (H) 

Use this freedom to impose the harmonic gauge condition h^ a > a = 0. Since 
the gauge remains unchanged for any transformation such that = 0, we 
can impose further conditions. We take one of these to be h£ = 0, to get 

W Q = D K» = o (12) 

That is, the metric variations are governed by a standard wave equation. 
Finally, use ft-oo = hjo = h j = to get hjk = hJ^F as before. The set of 
coordinates we have introduced is known as TT-gauge. 

Before we move on to discuss the modelling of various gravitational-wave 
sources, it is worth elucidating an issue that caused serious debate until the 
late 1960s. We need to demonstrate that gravitational waves carry energy. 
This is a tricky problem because, as we have already pointed out, one can 
only deduce the presence of a wave from the relative effect on two (or more) 
test particles. This means that one cannot localize the wave to individual 
points in space, and hence cannot directly "measure" its energy. In order 
to construct a meaningful energy expression we need to average over one 
(or more) wavelengths. Defining perturbations with respect to an averaged 
spacetime metric, i.e. using 

g$ = — 9,** = gfJ + V ( 13 ) 

and expanding in powers of h (which is presumed small), we have the 
(schematic) Einstein equations, 

<GW = G% + (G$)+ (G$) = (14) 



^o{h~) o(Taj 

Here, the term that is linear in h will vanish if we average over a wavelength. 
This means that we can deduce an expression for the stress energy-tensor for 
gravitational waves: 

^ = -ftr(f)^< W (15) 



Working out the algebra one can show that, in TT-gauge we will have 
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In particular, the energy propagating in the z-direction then follows from 

^ W - -^o W = -± (hi + hi)^E= -3£W = ^-(hl + hi) 

(17) 

where the frequency of the wave is lo. Finally assuming that h + ~ /i x ~ 
/isina>(t — z), and integrating over a sphere with radius r, we get 

E=^h>^\h\* = £LE (18) 
As we will now demonstrate, this is a very useful relation. 



2.2 Estimating the gravitational-wave amplitude 

We can use the expression (18) to infer the gravitational-wave strain asso- 
ciated with typical gravitational-wave sources. Let us characterise a given 
event by a timescale r and assume that the signal is monochromatic (with 
frequency /). Then we can use dE/dt w E/t and h w 2wfh to deduce that 



F, \ ' f t \ -i/2 / f 



1/2 

10 : >\[ .c 2 ) (l ms) \ i kHz / \ lo Mpc 

(19) 

If the signal analysis is based on matched filtering, the effective amplitude 
improves roughly as the square root of the number of observed cycles n. 
Using iiRi/rwc get 

h < - 5 x 1(r22 ( iF^) (r4) (15^) _1 ( 2 °) 

This is a crucial expression. We see that the "detector sensitivity" essentially 
depends only on the radiated energy, the characteristic frequency and the 
distance to the source. That is, in order to obtain a rough estimate of the 
relevance of a given gravitational-wave source at a given distance we only 
need to estimate the frequency and the radiated energy. Alternatively, if we 
know the energy released can work out the distance at which these sources 
can be detected. 

It is quite easy to obtain a rough idea of the frequencies involved. The 
dynamical frequency of any self-bound system with mass M and radius R is 



Thus, the natural frequency of a (non-rotating) black hole should be 
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(22) 



Medium sized black holes, with masses in the range 10 — IOOMq, will be 
prime sources for ground-based interferometers, while black holes with masses 
<~ 10 6 M Q should radiate in the LISA bandwidth. Meanwhile, neutron stars, 
with a mass of 1.4M compressed inside a radius of 10 km, will radiate at 



This means that one would expect neutron physics to be in the range for 
ground based detectors. In fact, given the likely need to detect signals with 
frequencies above 1 kHz, neutron star signals provide a strong motivation for 
the development of (third generation) high-frequency detectors. 

Given that the weak signals are going to be buried in detector noise, we 
need to obtain as accurate theoretical models as possible. The rough order 
of magnitude estimates we just derived will certainly not be sufficient, even 
though they provide an indication as to whether it is worth spending the 
time and effort required to build a detailed model. Such source models are 
typically obtained using either 

approximate perturbation techniques, eg. expansions in small perturba- 
tions away from a known solution to the Einstein equations, the archetypal 
case being black-hole and neutron star oscillations. 

post-Newtonian approximations, essentially an expansion in the ratio be- 
tween a characteristic velocity of the system and the speed of light, most 
often used to model the inspiral phase of a compact binary system, 
numerical relativity, where the Einstein equations are formulated as an 
initial- value problem and solved on the computer. This is the only way 
to make progress in situations where the full nonlinearitics of the theory 
must be included, eg. in the merger of black holes and neutron stars or a 
supernova core collapse. 

Here we will only describe the first step beyond a Newtonian description, 
where the gravitational radiation is described by the so-called quadrupole 
formula. For a source with weak internal gravity we have (in TT-gauge) 



We can solve this equation using the standard retarded Green's function to 
get 



Matching of the near-zone solution to an outgoing wave solution far away 
leads to the expression 



/ NS « 2 kHz 



(23) 



Uhp, = -167T1) 



(24) 




(25) 




(26) 



Gravitational-wave astronomy: the high-frequency window 9 

where 

T jk = J p (^XjXk - ^r 2 Sjk^ d 3 x (27) 

is the reduced quadrupole moment of the source. Consider a system of mass 
M with typical internal velocity v. Then we see that 

which shows (no surprise!) that in order to generate strong gravitational 
waves we need large masses moving at high speeds. 

From the formulas we derived earlier, we find that energy is radiated at 
a rate 

<»> 

The radiated angular momentum follows from the (usually) weaker current 
multipole radiation, which is governed by a similar expression. Let us now 



z 




Fig. 3. A schematic illustration of a compact binary system. 



apply the above results to the potentially most important gravitational-wave 
source, a compact binary system. Gravitational waves are emitted as the 
stars (or black holes) orbit each other and as a result the binary separation 
decreases. Consider a binary system with individual masses M\ and M 2 and 
separation 2R. Introduce the total and reduced masses 

Mi Mo 

M = M x +M 2 and fi = 1 — ► M Y a x = M 2 a 2 = (30) 

and work in the coordinate system illustrated in figure 3. Working out the 
required (time- varying) components of the quadrupole moment, we have 

1 XX = -T V y = cos 2<j) and T xy = 1 yx = sin 2</> (31) 



and wc find that 
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Next, determining the orbital rotation frequency from Kepler's law fP = 
M/i? 3 , and introducing the so-called "chirp mass" M. = [j, 3 / 5 M 2 / 5 we have 
the final result 

Moreover, from 
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Fig. 4. Estimated signal strengths for various inspiralling binaries relevant for 
ground- and space-based detectors. 



TT 

XX 



cos2f2(t - r) 



we can estimate the effective amplitude of the binary signal: 
h r ' 



(34) 



(35) 



This shows that, even though the actual signal gets stronger, its detectability 
decreases as the orbit shrinks. Figure 4 compares our estimated gravitational- 
wave strain to the predicted noise-curves for various gravitational-wave de- 
tectors. Essentially, one would expect that 

Advanced LIGO may observe several binary systems per year. 
The space-based LISA detector may suffer an "embarrassment of riches" , 
with a large number of known galactic binaries leading to detectable sig- 
nals (most likely generating to a "binary noise" which will be difficult to 
filter out). 
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To predict the rate at which the binary orbit shrinks as a result of 
gravitational- wave emission we need to estimate to total energy of the system: 

M lV \ M 2 v 2 2 GM X M 2 GfiM M 5/3 f2 2/3 

E = — + — — = -^m=- G 2 (36) 

From this we find that the period P of the system changes as 

For the binary pulsar 1913+16, the predicted change in orbital period agrees 
with the theoretical prediction to within 1%. This indirect proof that gravi- 
tational waves exist led to Hulse and Taylor being awarded the 1993 Nobel 
prize in physics. 

Finally, note that that the chirp-mass M. plays a key role in our various 
expressions. From 

t^ M ^ n ^ and h ~ M5/3nV3 (38) 
P r 

we see that this is the only combination of the two masses that can be infered 
from the signal (at this level of approximation: higher order post-Newtonian 
corrections depend on the individual masses, the spins etcetera). Suppose that 
we detect both the change in the amplitude h and the shift in the frequency. 
Then one can infer both the chirp mass and the distance to the source, and in 
effect coalescing binaries are standard candles which can be used to constrain 
cosmological parameters. 

One often classifies gravitational-wave sources by the nature of the waves. 
This is convenient because the different classes require different approaches 
to the data-analysis problem; 

Chirps. As a binary system radiates gravitational waves and loses energy 
the two constituents spiral closer together. As the separation decreases the 
gravitational-wave amplitude increases, leading to a characteristic "chirp" 
signal. 

Bursts. Many scenarios lead to burst-like gravitational waves. A typical 
example would be black-hole oscillations excited during binary merger. 
Periodic. Systems where the gravitational-wave backreaction leads to 
a slow evolution (compared to the observation time) may radiate per- 
sistent waves with a virtually constant frequency. This would be the 
gravitational- wave analogue of the radio pulsars. 

Stochastic. A stochastic (non-thermal) background of gravitational waves 
is expected to have been generated following the Big Bang. One may also 
have to deal with stochastic gravitational-wave signals when the sources 
are too abundant for us to distinguish them as individuals. 
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3 High-frequency gravitational wave sources 

Having introduced the key concepts required for a discussion of gravitational- 
wave physics we will now focus our attention on sources that radiate above a 
hundred Hz or so, i.e. which may at least in principle be detectable from the 
ground. As we will see there are strong motivations for constructing detectors 
which are sensitive up to (ideally) several kHz. 

3.1 Radiation from binary systems 

It is easy to understand why binary systems are considered the "best" sources 
of gravitational waves. They emit copious amounts or gravitational radiation, 
and for a given system we know quite accurately the amplitude and frequency 
of the gravitational waves in terms of the masses of the two bodies and their 
separation (see section 2.2). 

The gravitational-wave signal from inspiraling binaries is approximately 
sinusoidal, see equation (34), with a frequency which is twice the orbital 
frequency of the binary. As the binary system evolves the orbit shrinks and 
the frequency increases in the characteristic chirp. Eventually, depending 
on the masses of the binaries, the frequency of the emitted gravitational 
waves will enter the bandwidth of the detector at the low-frequency end 
and will evolve quite fast towards higher frequencies. A system consisting 
of two neutron stars will be detectable by LIGO when the frequency of the 
gravitational waves is ~10Hz until the final coalescence around 1 kHz. This 
process will last for about 15 min and the total number of observed cycles will 
be of the order of 10 4 , which leads to an enhancement of the dctectability 
by a factor 100 (remember h c <~ y/nh). Binary neutron star systems and 
binary black hole systems with masses of the order of 50Mq are the primary 
sources for LIGO. Given the anticipated sensitivity of LIGO, binary black 
hole systems are the most promising sources and could be detected as far as 
200 Mpc away. For the present estimated sensitivity of LIGO the event rate 
is probably a few per year, but future improvements of detector sensitivity 
(the LIGO II phase) could lead to the detection of at least one event per 
month. Supermassive black hole systems of a few million solar masses are 
the primary source for LISA. These binary systems are rare, but due to the 
huge amount of energy released, they should be detectable from as far away 
as the boundaries of the observable universe. Finally, the recent discovery of 
the highly relativistic binary pulsar J0737-3039 [4] enchanced considerably 
the expected coalescence event rate of NS-NS binaries [5] . The event rate for 
initial LIGO is in the best case 0.2 per year while advanced LIGO might be 
able to detect 20-1000 events per year. 

3.2 Gravitational collapse 

One of the most spectacular events in the Universe is the supernova (SN) 
collapse to create a neutron star (NS) or a black hole (BH). Core collapse is 
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a very complicated event and a proper study demands a deep understanding 
of neutrino emission, amplification of the magnetic fields, angular momen- 
tum distribution, pulsar kicks, etc. There are many viable models for each of 
the above issues but it is still not possible to combine all of them together 
into a consistent explanation. Gravitational waves emanating from the very 
first moments of the core collapse might shed light on all the above problems 
and help us understand the details of this dramatic event. Gravitational col- 
lapse compresses matter to nuclear densities, and is responsible for the core 
bounce and the shock generation. The event proceeds extremely fast, lasting 
less than a second, and the dense fluid undergoes motions with relativistic 
speeds (v/c ~ 0.2 — 0.4). Even small deviations from spherical symmetry 
during this phase can generate significant amounts of gravitational waves. 
However, the size of these asymmetries is not known. From observations in 
the electromagnetic spectrum we know that stars more massive than ~ 8M© 
end their evolution in core collapse and that ~ 90% of them are stars with 
masses ~ 8 — 20M Q . During the collapse most of the material is ejected and 
if the progenitor star has a mass M < 20M© it leaves behind a neutron star. 
If M > 20M Q more than 10% falls back and pushes the proto-neutron-star 
(PNS) above the maximum NS mass leading to the formation of a black hole 
(type II collapsars). Finally, if the progenitor star has a mass M > 40M© 
no supernova takes place. Instead, the star collapses directly to a BH (type I 
collapsars). 

A significant amount of the ejected material can fall back, subsequently 
spinning up and reheating the nascent NS. Instabilities can be excited during 
such a process. If a BH was formed, its quasi-normal modes (QNM) can be 
excited for as long as the process lasts. "Collapsars" accrete material during 
the very first few seconds, at rates ~ 1 — 2M Q /sec. Later the accretion rate is 
reduced by an order of magnitude but still material is accreted for a few tenths 
of seconds. Typical frequencies of the emitted gravitational waves are in the 
range l-3kHz for <~ 3 — 10M© BHs. If the disk around the central object has 
a mass ~ 1M© self-gravity becomes important and gravitational instabilities 
(spiral arms, bars) might develop and radiate gravitational waves. Toroidal 
configurations can be also formed around the collapsed object. Their instabil- 
ities and oscillations might be an interesting source of gravitational waves [7]. 
There is also the possibility that the collapsed material might fragment into 
clumps, which orbit for some cycles like a binary system (fragmentation in- 
stability [8]). 

The supernova event rate is 1-2 per century per galaxy [6] and about 5- 
40% of them produce BHs through the fallback material [9] . Conservation of 
angular momentum suggests that the final objects should rotate close to the 
mass shedding limit, but this is still an open question, since there is limited 
knowledge of the initial rotation rate of the final compact object. Pulsar 
statistics suggest that the initial periods are probably considerably shorter 
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than 20 ms. This strong increase of rotation during the collapse has been 
observed in many numerical simulations (see e.g. [10, 11]). 

Core collapse as a potential source of gravitational waves has been stud- 
ied for more than three decades (some of the most recent calculations can 
be found in [12, 13, 14, 15, 11, 16, 17, 18, 19]). All these numerical calcula- 
tions show that signals from Galactic supernova (d ~ lOkpc) are detectable 
even with the initial LIGO/EGO sensitivity at frequencies <lkHz. Advanced 
interferometers can detect signals from distances of 1 Mpc but it will be 
difficult with the designed broadband sensitivity to resolve signals from the 
Virgo cluster (<~15Mpc). The typical gravitational wave amplitude from the 
2D numerical simulations [11, 16] for an observer located on the equatorial 
plane of the source is 



where e ~ 1 is the normalized gravitational wave amplitude. The total energy 
radiated in gravitational waves during the collapse is < 10~ 6 — 10 -8 M Q c 2 . 
However, these numerical estimates are not yet conclusive, as important as- 
pects such as 3D hydrodynamics combined with proper spacetime evolution 
have been neglected. The influence of the magnetic fields have been ignored 
in most calculations. The proper treatment of these issues might not change 
the above estimates by orders of magnitude but it will provide a conclusive 
answer. There are also issues that need to be understood such as the pulsar 
kicks (velocities even higher than 1000 km/s) which suggest that in a frac- 
tion of newly-born NSs (and BHs) the process may be strongly asymmetric 
[20, 21, 22, 23, 24]. Also, the polarization of the light spectra in supernovae 
indicates significant asymmetries [25]. Better treatment of the microphysics 
and construction of accurate progenitor models for the angular momentum 
distributions are needed. All these issues are under investigation by many 
groups. 

Accretion Incuded Collapse (AIC) is also a possible source of high fre- 
quency gravitational waves. AIC takes place when a white dwarf (WD) ex- 
ceeds the Chandrasekhar limit due to accretion of material and begins to 
collapse. The cooling via neutrino emission does not reduce the heating sig- 
nificantly and the collapsing WDs reach appropriate temperatures for ignition 
of nuclear burning (Type la supernova). Estimates suggest that about O.1M 
material is ejected. Since the WD is pushed over the Chandrasekhar limit due 
to accretion, it will rotate fast enough to allow various types of instabilities 
[26]. The galactic rate of accretion induced collapse is about 10~ 5 /yr which 
means that AIC are about 1000 times rarer than core collapse SN. 

3.3 Rotational instabilities 

Newly born neutron stars are expected to rotate rapidly enough to be sub- 
ject to rotation induced instabilities. These instabilities arise from non- 




(39) 
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axisymmetric perturbations having angular dependence e lmv . Early New- 
tonian estimates have shown that a dynamical bar-mode (m — 2) instability 
is excited if the ratio (3 — T/W of the rotational kinetic energy T to the 
gravitational binding energy W is larger than (3d yn = 0.27. The instability 
develops on a dynamical time scale (the time that a sound wave needs to 
travel across the star) which is about one rotation period, and may last from 
1 to 100 rotations depending on the degree of differential rotation in the PNS. 
Another class of instabilities are those driven by dissipative effects such as 
fluid viscosity or gravitational radiation. Their growth time is much longer 
(many rotational periods) but they can be excited for significantly lower ro- 
tational rates, f3 > 0.14 in the case of the fundamental modes of oscillation 
of the star. 

3.4 Bar-mode instability 

The dynamical bar-mode instability can be excited in a hot PNS, a few 
milliseconds after the core-bounce, given a sufficiently large (3. It might also 
be excited a few tenths of seconds later, when the NS cools enough due to 
neutrino emission and contracts still further ((3 ~ 1/-R). The amplitude of the 
emitted gravitational waves can be estimated as h ~ MR 2 f} 2 /d, where M 
is the mass of the body, R its size, Q the rotational rate and d the distance 
from Earth. This leads to an estimate of the gravitational wave amplitude 

™ 23 (^(s4) 2 (^)^ <«» 

where e measures the ellipticity of the bar. Note that the gravitational wave 
frequency / is twice the rotational frequency £1. Such a signal is detectable 
only from sources in our galaxy or the nearby ones (our Local Group). If 
the sensitivity of the detectors is improved in the kHz region, then signals 
from the Virgo cluster may be detectable. If the bar persists for many (~ 10- 
100) rotation periods, then even signals from distances considerably larger 
than the Virgo cluster could be detectable, cf. Figure 5. The event rate is of 
the same order as the SN rate (a few events per century per galaxy): this 
means that given the appropriate sensitivity at frequencies between l-3kHz 
we might be able to observe a few events per year. The bar-mode instability 
may also be excited during the merger of NS-NS, BH-NS, BH-WD and even 
in type II collapsars (see discussion in [27]). 

In general, the above estimates rely on Newtonian hydrodynamics calcu- 
lations; GR enhances the onset of the instability slightly, /3dyn ~ 0.24 [28] and 
/?dyn may be even lower for large values of the compactness (larger M/R). 
The bar-mode instability may be excited for significantly smaller (3 if centrifu- 
gal forces produce a peak in the density off the sources rotational center [29]. 
Rotating stars with a high degree of differential rotation are also dynamically 
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unstable for significantly lower /? dyn > 0.01 [30, 31]. In this scenario the un- 
stable neutron star settles down to a non-axisymmctric quasistationary state 
which is a strong emitter of quasi-periodic gravitational waves 

h „ , v in -22 ( R** \ ( f f 100M P c \ M 1 ' 2 UU 

/lcff " 3xl ° (,30k^yl ^SOOH^J {—d-) M ™- (41) 

The bar-mode instability of differentially rotating neutron stars could be an 
excellent source of gravitational waves provided that the dissipation of non- 
axisymmetric perturbations by viscosity and magnetic fields is negligible. 
That this is the case is far from clear. Magnetic fields might actually enforce 
the uniform rotation of the star on a dynamical timcscalc and the persistent 
non-axisymmetric structure might not have time to develop at all. 

Numerical simulations have shown that the m = 1 one-armed spiral mode 
might become dynamically unstable for considerably lower rotational rates 
[29, 32]. This m = 1 instability depends critically on the softness of the 
equation of state (EoS) and the degree of differential rotation. 

3.5 CFS instability, f- and r-modes 

After the initial bounce, neutron stars may maintain a considerable amount 
of deformation. They settle down to an axisymmetic configuration mainly 
due to emission of gravitational waves, viscosity and magnetic fields. Dur- 
ing this phase QNMs are excited. Technically speaking, an oscillating non- 
rotating star has equal values ±|cr| (the frequency of a mode) for the forward 
and backward propagating modes (corresponding to m = ±|m|). Rotation 
changes the mode frequency by an amount 5a ~ mfl and both the prograde 
and retrograde modes will be dragged forward by the stellar rotation. If the 
star spins sufficiently fast, the originally retrograde mode will appear to be 
moving forwards in the inertial frame (according to an observer at infinity), 
but still backwards in the rotating frame (for an observer rotating with the 
star). Thus, an inertial observer sees gravitational waves with positive an- 
gular momentum emitted by the retrograde mode, but since the perturbed 
fluid rotates slower than it would in absence of the perturbation, the angu- 
lar momentum of the mode itself is negative. The emission of gravitational 
waves consequently makes the angular momentum of the mode increasingly 
negative leading to an instability. From the above, one can easily conclude 
that a mode will be unstable if it is retrograde in the rotating frame and 
prograde for a distant observer measuring a mode frequency a — mfl i.e. the 
criterion will be a (a — mf2) < 0. 

This class of frame- dragging instabilities is usually referred to as Chandra- 
sckhar-Friedman-Schutz (CFS) [34, 35] instabilities. For the high frequency 
(/ and p) modes the instability is possible only for large values of Q or for 
quite large m. In general, for every mode there will always be a specific value 
of fi for which the mode will become unstable, although only modes with 
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|m| < 5 have an astrophysically significant growth time. The CFS mechanism 
is not only active for fluid modes but also for the spacetime or the so-called 
w- modes [33]. It is easy to see that the CFS mechanism is not unique to 
gravitational radiation: any radiative mechanism will have the same effect. 

In GR, the f-mode (I = m = 2) becomes unstable for (3 ~ 0.06 — 0.08 
[36]. If the star has significant differential rotation the instability is excited 
for somewhat higher values of (3 (see e.g. [37, 38]). The /-mode instabil- 
ity is an excellent source of gravitational waves. After the brief dynamical 
phase, the PNS becomes unstable and the instability deforms the star into 
a non-axisymmetric configuration via the I = 2 bar mode. Since the star 
loses angular momentum it spins down, and the gravitational wave frequency 
sweeps from 1 kHz down to about 100 Hz [39]. If properly modelled such a 
signal can be detected from a distance of 100 Mpc (if the mode grows to a 
large nonlinear amplitude). 

Rotation not only shifts the frequencies of the various modes; it also gives 
rise to the Coriolis force, and an associated new family of rotational or inertial 
modes. Inertial modes are primarily velocity perturbations. Of special interest 
is the quadrupole inertial mode (r-mode) with I = m = 2. The frequency of 
the r-mode in the rotating frame of reference is a = 2f2/3. Using the CFS 
criterion for stability we can easily show that the r-mode is unstable for 
any rotation rate of the star. For temperatures between 10 7 — 10 9 K and 
rotation rates larger than 5-10% of the Kepler limit, the growth time of the 
unstable mode is much shorter than the damping times due to bulk and 
shear viscosity. The mode grows until it saturates due to non-linear effects 
[40, 41, 42]. The strength of the emitted gravitational waves depends on the 
saturation amplitude a. Mode coupling might not allow the growth of the 
instability to amplitudes larger than a s=a 10~ 2 — 10~ 3 [45]. The existence of 
a crust [43, 44] or of hyperons in the core [47] and strong magnetic fields [46] , 
affect the efficiency of the instability (for extended reviews see [48, 49]). For 
newly-born neutron stars the amplitude of gravitational waves might not be 
such that the signals will be detectable only from the local group of galaxies 
{d < IMpc) 



see Figure 5. 

If the compact object is a strange star, then the r-mode instability will 
not reach high amplitudes (a <~ 10~ 3 — 10~ 4 ) but it will persist for a few 
hundred years and in this case there might be up to ten unstable stars per 
galaxy radiating gravitational waves at any time [54] . Integrating data for a 
few weeks can then lead to an effective amplitude /i c ff ~ 10~ 21 for galactic 
signals at frequencies ~ 700 — 1000Hz. The frequency of the signal changes 
only slightly on a timescale of a few months, so the radiation is practically 
monochromatic. 




(42) 
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Old accreting neutron stars, radiating gravitational waves due to the r- 
mode instability, at frequencies 400-700Hz, are probably a better source [50, 
51, 52, 53]. Still, the efficiency and the actual duration of the process depends 
on the saturation amplitude a. If the accreting compact object is a strange 
star or has a hyperon core then it might be a persistent source which radiates 
gravitational waves for as long as accretion lasts [54, 55] . 




gravitational-wave frequency (Hz) 



Fig. 5. The estimated strength of gravitational waves from the dynamical bar- 
mode instability and the CFS instability of the f- and r-modes. The estimates are 
compared to the predicted noise of the various interferometers and also the possible 
noise curve for a dual cylinder detector. 



3.6 Oscillations of black holes and neutron stars 

Black-hole ringing. The merger of two neutron stars or black holes or the 
collapse of a supermassive star (collapsar of type I or II) will produces a 
black hole. The newly formed black hole will ring, emitting a characteristic 
signal until it settles down to the stationary Kerr state. This characteristic 
signal, the so-called quasi-normal mode oscillation, will be a unique probe 
of the black hole's existence. Although the ringing phase does not last very 
long (a few tenths of a ms), the ringing due to the excitation by the fallback 
material might last for secs[56, 57]. The frequency and damping time of the 
black-hole ringing for the I = m = 2 oscillation mode can be estimated via 
the relations [58] 
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cr « 3.2kHz Mfp 1 1 - 0.63(1 - a/M) 3/10 (43) 

Q = irar w 2(1- a) ~ 9/20 (44) 

These relations together with similar ones either for the 2nd QNM or the 
I = 2, m = 0, ±1 can uniquely determine the mass M and angular momentum 
a of the BH if the frequency and the damping time of the signal have been 
accurately extracted [59, 60]. The amplitude of the ring-down waves depends 
on the BH's initial distortion. If the excitation of the BH is due to infalling 
material then the energy is roughly AE > e/j,c 2 (fi/M) where e > 0.01 [61]. 
This leads to an effective gravitational wave amplitude 

"(o*)(^)(£) 

This approximate result has been verified by more detailed full non-linear 
simulations [62] 

Neutron star ringing. If the collapse leaves behind a compact star, various 
types of oscillation modes might be excited which can help us estimate pa- 
rameters of the star such as radius, mass, rotation rate and EoS[63, 64, 65]. 
This gravitational wave asteroseismology is a unique way to find the radius 
and the EoS of compact stars. One can derive approximate formulas in or- 
der to connect the observable frequencies and damping times of the various 
stellar modes to the stellar parameters. For example, for the fundamental 
oscillation (I = 2) mode (/-mode) of non-rotating stars we get [63] 

cr(kHz) w 0.8 + l.&Ml'^R^' 2 + 5 im fi (46) 
T^sccs- 1 ) w Mf A R^ (22.9 - U.1M 1A R^) + 8 2 mfi (47) 

where Q is the normalized rotation frequency of the star, and 8\ and 82 
are constants estimated by sampling data from various EoS. The typical 
frequencies of the NS modes are higher than 1 kHz. On the other hand, 2D 
simulations of rotating core-collapse have shown that if a rapidly rotating 
NS is created, then the dominant mode is the quasi-radial mode ( u l — 0" ) , 
radiating through its I = 2 piece at frequencies <~ 800Hz-lkHz [11]. Since 
each type of mode is sensitive to the physical conditions where the amplitude 
of the mode eigenfunction is greatest, the more information we get from the 
various classes of modes the better we will understand the details of the star. 

Concluding, we should mention that the tidal disruption of a NS by a BH 
[66] or the merging of two NSs [67] may give valuable information for the 
radius and the EoS if we can detect the signal at frequencies higher than 1 
kHz. 
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